Phase competition in the double-exchange model on the frustrated pyrochlore lattice 
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Competition between the ferromagnetic double- exchange interaction and the super-exchange an- 
tiferromagnetic interaction is theoretically studied in the presence of geometrical frustration. As 
increasing the super-exchange interaction, the ferromagnetic metal becomes unstable, and is taken 
over by a cooperative paramagnetic metal, in sharp contrast with a discontinuous transition to the 
antiferromagnetic insulator in the absence of frustration. In the critical region, the system exhibits 
a peculiar temperature-independent behavior with highly incoherent transport, suggesting a large 
residual entropy at low temperatures. We discuss the relevance of the results to the pressure-induced 
behaviors in Mo pyrochlore oxides [S. Iguchi et al, Phys. Rev. Lett. 102, 136407 (2009)]. 

PACS numbers: 71.30.+h, 71.27.-fa, 71.20.Be, 71.10.Fd 



Interplay between spin and charge degrees of freedom is 
a central issue in the field of strongly-correlated electron 
systems. A main goal is to clarify how magnetic and 
charge fluctuations are related with each other and how 
various phase transitions such as magnetic ordering and 
metal-insulator transition emerges from the interplay. 

One of the fundamental models describing the inter- 
play explicitly is the double-exchange (DE) model [l|, 
which has been extensively studied for understanding 
the physics of colossal magnetoresistance (CMR) in per- 
ovskite manganese oxides [H. It has been clarified that 
the model exhibits a first-order phase transition with a 
bicritical behavior between a ferromagnetic metal (FM) 
and an antiferromagnetic (AF) insulator \§, The 
FM is stabilized by the DE interaction which is an ef- 
fective ferromagnetic interaction originating from the 
strong coupling between itinerant electrons and local- 
ized moments, while the AF insulator is stabilized by 
the super-exchange (SE) AF interaction between local- 
ized moments. It was also discussed that the FM state 
competes with a charge- and/or orbital-ordered insula- 
tor in extended DE models including the electron-lattice 
coupling and the orbital degeneracy [1, 0] . 

In general, magnetism is strongly infiuenced by the ge- 
ometry of underlying lattice structure. In the absence 
of geometrical frustration, AF interaction leads to a sim- 
ple Neel-type AF ordering, whereas it becomes unstable 
against the frustration, resulting in nontrivial phenom- 
ena such as a complicated ordering, a glassy state, and 
a spin-liquid state 0- The frustrated magnetism in lo- 
calized spin systems is an interesting long-standing issue, 
but it will be more intriguing to clarify the effect of geo- 
metrical frustration in spin-charge coupled systems. The 
interplay between the frustrated magnetism and the elec- 
tronic transport is expected to yield characteristic spin- 
charge entangled phenomena. 

There are many compounds in which the spin and 
charge degrees of freedom couple with each other un- 
der frustration. Among them. Mo pyrochlore oxides 
i?2Mo2 07 are interesting because they exhibit phase 
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FIG. 1: (Color online) Phase diagram of the pyrochlore 
double-exchange model at quarter filling. FM, PS, PM, and 
CPM represent ferromagnetic metal, phase separation, para- 
magnetic metal, and cooperative paramagnetic metal, respec- 
tively. Closed symbols represent the phase transition bound- 
aries, while open ones indicate the magnitude of Curie- Weiss 
temperature. The lines are guides for the eye. The inset shows 
the cubic unit cell of the pyrochlore lattice. The dotted bond 
shows an example of the ID chain for the calculations of the 
spin correlation in Fig. [31 



competition between FM and spin-dass insulator (SGI) 
by changing the rare earth R [1, In these com- 
pounds. Mo cations constitute the severely-frustrated 
pyrochlore lattice, which consists of the 3D network of 
corner-sharing tetrahedra, as shown in the inset of Fig. [T] 
It was pointed out by the first-principle band calculation 
that FM is stabilized by the DE interaction, while the 
transition to SGI is triggered by the Coulomb repulsion 
between Mo 4d electrons and the resulting Mott localiza- 
tion [m. 

Recently, several new features were revealed in the Mo 
pyrochlore oxides under external pressure [Til . [T2j . The 
pressure induces a transition from FM to a peculiar para- 
magnetic metallic (PM) state, in which the resistivity is 
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highly incoherent and almost teniperature(r) indepen- 
dent. In addition, a spin-glass metallic (SGM) phase is 
found in between. There is no metal-insulator transition, 
which is clearly distinguished from the situation in the R- 
site substitution at ambient pressure. It is highly desired 
to elucidate relevant interactions to the pressure-induced 
behaviors. 

Motivated by both the theoretical and the experimen- 
tal problems, in the present study, we investigate the 
phase competition in the DE model defined on the 3D 
frustrated pyrochlore lattice. We employ the unbiased 
Monte Carlo simulation and illuminate distinctive as- 
pects of the frustrated spin-charge coupled system. 

The Hamiltonian of the DE model is given by 



{ij)a i 



Jaf ^ Si 



(1) 

where t is the transfer integral for nearest-neighbor (n.n.) 
sites (ij) on the pyrochlore lattice in the inset of Fig. [H 
Jh represents the Hund's-rule coupling between the itin- 
erant electron spin and the localized spin Si, and Jaf 
is the SE interaction between the neighboring localized 
spins. Some possible symmetry-broken states were ex- 
amined for this model by mean-field approximation jisj . 
For simplicity, hereafter, we consider the limit oi Ju ~^ oo 
and treat the localized spins Si as classical vectors with 
a normalized length | Si | = 1 . We set an energy unit 
as i = 1. In the following, we investigate the properties 
of the mode l Jp at quarter filling (0.5 electron per site 
on average) |15| . 

We employ the Monte Carlo simulation to study the 
thermodynamic properties of the model The algo- 
rithm is a standard one in which configurations of clas- 
sical localized spins are sampled by Monte Carlo proce- 
dure; the Monte Carlo weight is calculated by the exact 
diagonalization of the fermion Hamiltonian matrix for a 
given spin configuration Tg']. The system sizes Ng are 16, 
32, 64, and 128 sites which correspond to Ix 1x1, Ix 1x2, 
2 X 2x1, and 2x2x2 cubic unit cells. We take an aver- 
age over the twisted boundary conditions [13, [3 , which 
enable systematic analysis within the small size clusters 
down to r ~ 0.01. The details and the efficiency of this 
technique are found in Ref. [19i] . 

Figure [1] shows the phase diagram of the model H]) at 
quarter filling. The ferromagnetic transition temperature 
Tc decreases gradually as increasing Jaf , and approaches 
zero at Jaf '^0.1. Here, Tq is estimated by a systematic 
analysis of the Binder parameter [19]. After the FM state 
fades out, we do not find any clear sign of phase transi- 
tion as will be discussed below. (We will mention about 
|8| and the CPM state later.) These behaviors are con- 
trastive to the results in the absence of frustration, where 
a bicritical behavior appears with a first-order transition 
to an AF insulating state [3l,l3|. In the present model, the 
frustration of the pyrochlore structure is severe enough 
to suppress the AF ordering and replace it with a para- 
magnetic state. 
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FIG. 2: (Color online) Temperature dependences of (a) the 
squared magnetic moment and (b) the inverse of the uni- 
form magnetic susceptibility, (c) Curie- Weiss temperature 
estimated from the fitting of x for various system sizes A^s- 
(d) Temperature dependence of the uniform magnetic suscep- 
tibility for large Jaf. In (a), (b), and (d), circles and crosses 
are the data for A^s = 128 and 64, respectively. 



Instead of the bicritical behavior in the unfrustrated 
cases, an electronic phase separation (PS) takes place at 
low T as shown in the phase diagram. In this region, the 
system is phase-separated into FM at a lower density and 
FM at a higher density, that is, quarter filling is no longer 
stable: PS is identified by a jump of the electron density 
as a function of the chemical potential. This metal-to- 
metal phase separation is a characteristic feature due to 
the frustration which suppresses AF ordering, however, 
we focus on the competition between FM and PM at 
quarter filling in the following, except for a brief comment 
related to experiments later. The details of PS will be 
discussed elsewhere. 

In Figs.[2Ua) and[2jb), we show T dependences of the 
squared moment = Si/A's)^) and the inverse of 
the uniform magnetic susceptibility, x^^- In ths small 
Jaf region, the ferromagnetic moment grows rapidly as 
decreasing T, and at the same time y shows a peak whose 
height grows with the system size These signal the 
ferromagnetic transitions. As Jaf increases, however, 
the growth of is suppressed gradually, corresponding 
to the collapse of FM in Fig. [T] In the larger- Jaf region, 
becomes zero and x is suppressed. 

As shown in Fig. ^h), x exhibits a Curie- Weiss like 
behavior at high T; x^^ cx T — O, where Q is the Curie- 
Weiss temperature. The estimate of 8, obtained by fit- 
ting the high-T part, is plotted in Fig. [^Jc). Q decreases 
from positive to negative continuously as increasing Jaf, 
and crosses zero at Jaf '^0.1. This indicates that the ef- 
fective magnetic interaction changes from ferromagnetic 
to antiferromagnetic continuously, and vanishes near the 
critical region at Jaf '^0.1 where Tq goes to zero. This 
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FIG. 3: (Color online) Spin correlation (Si -Sj) along the ID 
chains in the pyrochlore structure (see the inset of Fig. [T]) as a 
function of distance (measured in unit of the nearest- neighbor 
bond length) for (a) Jaf = 0.0, (b) 0.1, and (c) 0.2. The data 
are for iV^ = 128. 



suggests a cancellation between the DE ferromagnetic in- 
teraction and the SE AF interaction. 

Such cancellation is also observed in the spin correla- 
tion. Figure [3] shows the distance- dependence of the spin 
correlation (Si ■ Sj) for three typical values of Jaf = 0.0, 
0.1, and 0.2, corresponding to FM, critical, and PM re- 
gions, respectively. In the FM region [Fig.[31Ja)], the spin 
correlation grows rapidly as decreasing T, and converges 
to a positive value for large distance. In the critical re- 
gion at Jaf ~ 0.1, however, the spin correlation hardly 
develops as shown in Fig. [31[b). This indicates that the 
spins remain disordered due to the cancellation between 
the DE ferromagnetic interaction and the SE AF inter- 
action. 

In the PM region with large Jaf, AF correlation de- 
velops only for short range at low T [Fig. |3Ic)]. This is 
remarkable since Q increases its magnitude with Jaf, as 
shown in Fig.[5l^c), indicating a large effective AF interac- 
tion Jeff ^ 6/8. The suppressed AF correlation as well as 
no clear sign of phase transition is due to the strong frus- 
tration inherent to the pyrochlore lattice structure. This 
type of frustrated paramagnetic state below T ~ \Q\ is 
often called a cooperative paramagnet or a classical spin 
liquid [2^ [2l[ . We call the region as cooperative param- 
agnetic metal (CPM) and indicate it in Fig.[T]by plotting 
9|. Note that the susceptibility in the CPM region ex- 
hibits a characteristic behavior at low T; it deviates from 
the high-T Curie- Weiss behavior and shows a broad peak 
at T ~ Jeff, as plotted in Fig. [UJd). 

These drastic changes in the magnetic sector affect the 
electronic state through the strong spin-charge coupling 
in the model ([ij. In Fig. SI we present the density of 
states (DOS) D{lli) and the optical conductivity <j{lij) for 
the same set of Jaf in Fig. [31 DOS is plotted by setting 
the chemical potential to be zero, and (t{uj), a diagonal el- 
ement along the cubic axis, is calculated by the standard 
Kubo formula. In the FM region [Figs. [4][a) andlUJd)], 
both D{uj) and 17(0;) show a strong T dependence, cor- 
responding to the development of ferromagnetism. Note 
that the low-T results are similar to those for the nonin- 
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FIG. 4: (Color onhne) (a)-(c) Density of states and (d)-(f) 
optical conductivity. The data are for A's = 128. 



teracting model (two dispersive bands for —8 < w < —4 
and —4 < w < 0, and two degenerate flat bands at w = 0) 
because of less spin scattering from almost fully-aligned 
moments. (The spiky structure seen in the dispersive 
bands is due to the finite-size effect.) On the other hand, 
in the large- Jaf PM region [Figs, (tjc) and|li;f)[, T de- 
pendence is suppressed but still remains at low T, reflect- 
ing the growth of short-range AF correlations in Fig.[3Jc) 
[22| . In the critical region, neither D{ijj) nor <t{lij) exhibits 
any significant T dependence as exemplified in Figs.HJb) 
andjU^e), corresponding to the suppression of spin cor- 
relation down to the lowest T in Fig. |2Ib). Here the 
electrons are scattered by almost uncorrelated moments, 
resulting in the highly-incoherent conductivity shown in 
Fig. life). Thus our model exhibits a T-independent, in- 
coherent electronic state due to the competition between 
the DE ferromagnetic interaction and the SE AF inter- 
action. 

In this critical region, spin fiuctuations remain down 
to low T, leading to a large amount of residual en- 
tropy. In Fig. [Sja), we plot the internal energy per site, 
E = {H) /Ng, as a function of T. For small Jaf, E shows 
a rapid decrease as T decreases and an infiection point 
(indicated by arrows in the figure), corresponding to the 
ferromagnetic transition, but the decrease is gradually 
suppressed as Jaf increases. At Jaf ~ 0.1, the T de- 
pendence of E becomes minimum, whereas it is slowly 
recovered for larger Jaf without showing any anomaly 
associated with a phase transition. Hence E shows the 
smallest T dependence in the critical region. Here both 
the kinetic energy of electrons and the SE energy between 
localized moments become less T dependent, as shown in 
Figs.Hfb) and|5{c). These results suggest a large amount 
of entropy remaining at lower T, since the entropy is 
given by T integration of the specific heat dE/dT di- 
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FIG. 5: (Color online) Temperature dependence of (a) the 
total internal energy, (b) the kinetic energy, and (c) the SE 
energy. Circles and crosses represent the data for A's = 128 
and 64, respectively. In (a), arrows indicate inflection points 
which correspond to the transition to FM state (see text). 

vided by T. Although the entropy is not well defined 
in the model including classical vector spins, we believe 
that the almost T-independent behavior is robust even 
for more realistic models and that the residual entropy 
is widely observed in the frustrated DE-type models. If 
this PM state survives down to the lowest T, the sys- 
tem will be a non Fermi liquid [IS US] and show a large 
mass behavior. Another possibility is that some symme- 
try breaking eventually takes place at a lower T beyond 
our numerical limitation. 

Finally let us discuss the implication of our results to 
experiments. Mo pyrochlore oxides i?2Mo207 are found 
to show puzzling behaviors under external pressure in 
contrast to the 7?-site substitution at ambient pressure, 
that is, a collapse of FM and an emergence of the peculiar 
PM state with showing T-independent incoherent trans- 
port 12]. These are well reproduced in our results. In 
addition, the experimentally-observed SGM state may be 
accounted by our PS state: Coexisting phases potentially 
lead to domain formation in the presence of long-range 
interactions such as Coulomb repulsion, and the domains 
pinned by randomness can result in cluster-glass like be- 
haviors '25"|. The agreement suggests that (i) our simple 
DE model captures an essential physics of the pressure ef- 



fect on the complicated compounds and (ii) the geometri- 
cal frustration plays a major role in the pressure-induced 
phenomena where competition between DE and SE in- 
teractions are important. These are in sharp contrast 
to the previous study which focused on the relevance of 
Coulomb repulsion in the i?-site substitution at ambient 
pressure where the metal-insulator transition is observed 
[lO|. Such observations are important for understanding 
of the distinctive physics between the i?-site substitution 
and the external pressure. Further experiment on the 
structural change under pressure will help to understand 
the microscopic origin of the difference between the two 
effects. Our result, moreover, predicts a non-Fermi-liquid 
behavior with a large quasiparticle mass or some symme- 
try breaking at low T in the critical region. It is highly 
desired to examine experimentally the lower T behavior 
under pressure in more detail. 

To summarize, we have studied the phase competi- 
tion in the double-exchange model on the frustrated 
pyrochlore lattice. In sharp contrast with the unfrus- 
trated models studied for the CMR manganites which 
exhibit bicritical behaviors between ferromagnetic metal 
and antiferromagnetic insulator, we found that the fer- 
romagnetic metal gradually collapses against the super- 
exchange interaction, and finally a paramagnetic incoher- 
ent metal takes over. In the critical region, the localized 
moments are almost uncorrelated down to low tempera- 
tures, and the electronic transport becomes almost tem- 
perature independent and highly incoherent. The results 
are favorably compared with the pressure-induced phe- 
nomena in Mo pyrochlore oxides, and in addition, predict 
a large residual entropy. These interesting aspects orig- 
inate from the interplay between itinerant electrons and 
frustrated magnetism. Such interplay is widely seen in 
many frustrated materials and awaits future study. 
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